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Abstract 

Using Ito's calculus and the mass optimal transportation theory, we study 
the generalized Dyson Brownian motion (GDBM) and the associated McKean- 
Vlasov evolution equation with an external potential V. Under suitable condi- 
tion on V, we prove the existence and uniqueness of strong solution to SDE for 
GDBM. Standard argument shows that the family of the process of empirical 
measures L^{t) of GDBM is tight and every accumulative point of £jv(i) in the 
weak convergence topology is a weak solution of the associated McKean-Vlasov 
evolution equation, which can be realized as the gradient flow of the Voiculcscu 
free entropy on the Wasserstein space over R. Under the condition V" > — K 
for some constant K > 0, we prove that the McKean-Vlasov equation has a 
unique solution [i(t) and Ljy(t) converges weakly to fi(t) as N — > oo. For C 2 
convex potentials, we prove that /x(t) converges to the equilibrium measure 
fiv with respect to the W%- Wasserstein distance on .^(R) as t — >• oo. Un- 
der the uniform convexity or a modified uniform convexity condition on V, we 
prove that /j,(t) converges to /iy with respect to the ^-Wasserstein distance on 
^2(R) with an exponential rate as t — > oo. Finally, we discuss the double- well 
potentials and raise some conjectures. 

Key words and phrases: Generalized Dyson Brownian motion, McKean-Vlasov equation, 
Johansson's theorem, gradient flow, Voiculescu free entropy, Wasserstein distance. 

1 Introduction 

1.1 Background 

Let V : M — > R + be a real polynomial of even degree with positive leading coefficient, or 
more general a real analytic function. Consider the following probability measure on %m 
(the set of all N x N Hermitian matrices) 

dfJ, N (M) = — exp(- NTrV{M))dM, 
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where Zm is a normalization constant, and if we denote x±, 



. , xn the eigenvalues of M, 



N 



TiV{M)=Y J V{x t ). 

i=l 

By [7, 20] , the distribution of the eigenvalues of M has the following probability density 

1 ( N \ 

Pn{x) = — Ui<j\xi - xj\ 2 exp -Nj~]V(xi) , x 

Zn V ti J 



eR N . 



The probability distribution pn(x)<Ix has a statistical mechanical interpretation: it is 
the canonical Gibbs measure, at inverse temperature /3 = 2, of a system of N unit charges 
interacting through the logarithmic Coulomb potential and confined by an external potential 
NV. From the statistical mechanical point of view, it is natural to consider the logarithmic 
Coulomb gas at arbitrary values of the inverse temperature j3 > 0. More generally, let 
V : M — > R be a polynomial of even degree or an analytic function with the growth condition 

V(x) > (l + <5)log(x 2 + l), xeM, (1) 

the /3-invariant random matrix ensemble, or the so-called log-gas model, is defined as an 
interacting particle system with the following probability distribution 

1 p ( f3N N \ N 

Pp (dxi,. . -,dx N ) = —Hi^j\xi -Xj\* exp I — ^-^V^) J JJda;,, 

N \ i=l J i=l 

where /? > is a parameter, called the inverse of the temperature in statistical physics. When 
(3 = 1,2,4, one can realize the above distribution as the distribution of the eigenvalues of 
N x N random matrices. More precisely, for f3 — 1,2,4, if M is a N x N matrix with 
distribution 



dP»(M) = — exp ^-^-Try(M)J dM, 
where dM denotes the standard measure on Hn, i-e., 

dM = II dMu II dRe(M t3 )dlmM tj . 

l<i<N l<i<j<N 

In particular, for V(x) = y, we get Gaussian Symmetric matrices ensemble (GOE) for 
/3 = 1, the Gaussian Hermitian matrices ensemble (GUE) for (3 = 2, and for the Gaussian 
Symplectical matrices ensemble (GSE) for /? = 4. 

The distribution of interacting particles with general external potential V and the loga- 
rithmic Coulomb interaction has received a lot of attention in theoretic physics in connection 
with the so-called matrix models, see e.g. [9, 17]. We would like also to mention that, Kont- 
sevich [21] also used the complex partition function for the matrix model with V(x) = x 3 on 
Hermitian random matrices to prove Witten's conjecture in the intersection theory of the 
moduli space of curves. 
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Suppose that V : R — > R satisfies (1). In [20], Johansson proved the following result: 
There is a unique equilibrium measure /iy € <^(R) with compact support such that 



inf E v [/j] = Ey[/iv], 



and satisfies the Euler-Lagrange equation 

P.V. / = -V (a;), x e supply, 

Jr x ~ y £ 

where ^(R) is the set of all Bore probability measures on R, and £y is the following energy 
functional, the so-called Voiculescu free entropy functional 



5V(m) = - / / log | a; - y\d/j,(x)d(j,(y) + / V(x)dn(x). 
Jr Jr Jr 

n 

Moreover, as N tends to infinity, the expectation of the empirical measure Ljy = ^ y] 6\ i 

weakly converges to /iy, i.e., E [Ln] — > My- This recovers Wigner's famous semi-circle law 
[42] for GUE with V{x) = 

Let us consider the Gaussian Hermitian Ensemble (GUE). Let B t be the NxN Hcrmitian 
matrice with entries Bij(t), where Bij(t), 1 < i < j < N, arc i.i.d. complex Brownian 
motions with E[Bij(t)] = and E[|Sjj | 2 (i)] = t. Let Ai(i), . . . , Ajv(i) be the process of 
the eigenvalues of B t . In [15], Dyson proved that Xi(t) satisfies the following stochastic 
differential equations 

dXUt) = -^dWf + — dt. 

See also Mehta [28], Guionnet [18], Andersson-Guionnet-Zeitouni [2] and Tao [38]. 

The Dyson Brownian motion {Xi(t),i = 1, . . . ,7V} is an interacting iV-particle system 
with the logarithmic Coulomb interaction. It has been very useful in various branches of 
mathematics and physics, including statistical physics and the quantum chaotic systems. 
See the references mentioned in [28, 18, 2, 34]. Instead of using Hermitian matrix- valued 
Brownian motion, Chan [12], Rogers and Shi [36] proved that the eigenvalues of the NxN 
Hermitian Ornsetin-Uhlenbeck process satisfies the following stochastic differential equations 



dX 



N 



N 

They [12, 36] proved the tightness of the family of the empirical measure Lpj(t) = Aj(i) 

i=l 

for the above Dyson-Ornstein-Uhlenbeck Brownian motion (Ai(i), . . . , Ajv(i)) and proved 
that the limit of any weakly convergence subsequence of Ljv(i) is a weak solution of the 
following McKean-Vlasov equation 

4 / f{x)n t {dx) = \ ( ( dxf{x) - d » m K(dx)iH(dy)-± ( xf(x)to(dx), 
at J R 2j J R 2 x-y 2 J R 
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where / is a test function in C%(M. N ). In particular, the Hilbert-Stiejies transformation of 
n(t), defined by 



Gt(z) 



Ht{dx) 



satisfies the following nonlinear Burgers type equation 

l Gt (z) = (-G t (z) + z)^- z G t (z) + G t (z). 

Showing that the above complex Burgers equation has a unique solution, they derived that 
L^it) converges weakly to fi(t). Moreover, as t — > oo, Gt(z) converges to the Hilbert-Stiejies 
transform of the semi-circle law, one can therefore give a new proof of Wigner's theorem 
[42] for the convergence of the empirical measure of the eigenvalues of Gaussian Unitary 
Ensemble (GUE) to the semi-circle law. See also Guionnet [18] and Andersson-Guionnet- 
Zeitouni [2] for a nice presentation of this dynamic approach to Wigner's theorem using 
Dyson's Bownian motion. 

1.2 Motivation 

The purpose of this paper is to study the generalized Dyson Brownian motion and related 
McKean-Vlasov equation associated with the log-gas model with non-quadratic external 
potential. Thus, we are working on non Gaussian type /3-invariant ensembles with a general 
external potential V. To describe our motivation, let us first introduce the generalized Dyson 
Brownian motions (briefly, GDBM) as follows. 

Let (W 1 , . . . , W N ) be an ./V-dimensional Brownian motion defined on a probability 
space (fi,P) with a filtration & = {& t ,t > 0}. Let Ajv(O) = (A^(0), . . . , A#(0)) e A N , 
where 

A w = {(xi)i<i<Ar e K. N : xi < x 2 < ■ . . < x N }. 

By Theorem 1.2 below, for a wide class of potential function V with suitable condition, the 
following stochastic differential equations 



dx^t) - sj^wi + 1 g ^v-W) * " \ V ' [X ^ i = 1 >- N > (2) 

have a unique strong solution (Ajv(i))t>o with initial data Aat(0) such that Ajv(t) G Ajv 
for all t > 0. 

The process (Ajv(f))t>o defined by SDE (2) is called the generalized Dyson Brownian 
motion (GDBM) with potential V. The GDBM is an interacting particle system with 
Hamiltonian of the form 

1 1 N 

H(xi,...,x N ):=-— l °S\x l -x :i \ + -^2v(x l ). 

l<i^j<N i=l 



The infinitesimal generator of these interacting particles is given by 



dm 

dx k 



4 



where / e C 2 (R N ) and L N = ± £ S Xi e 5»(R). 

In particular, when V = 0, (Ajv(£))t>o is the Dyson Brownian motion introduced by F. 
Dyson [14, 15], and when V(x) = it is the Dyson-Ornstein-Uhlenbeck Brownian motion 
introduced by Chan [12] and Rogers and Shi [36]. In fact, as the same situation as what 

2 

has been done by Dyson for V = 0, and by Chan [12] and Rogers-Shi [36] for V(x) = \, 
we can prove that, in the case f3 = 1, 2,4, the generalized Dyson Brownian motion can be 
realized as the process of the eigenvalues of a matrix- valued diffusion process. See Remark 
1.6, Section 2.2 and Section 2.3 below. 

The process of the empirical measures of GDBM is defined by 

1 N 

Mi) = ^5>A. N (t)e^(iR), *e[o,oo). 

i=l 

In Section 3, we will prove that, under suitable condition on V, the family of the process of 
empirical measures Lpj(t) is tight, and any weak convergence limit of convergent subsequence 
of £jv(i), denoted by /i(t), is a weak solution to the following nonlinear McKean-Vlasov 
equation: for all / e Cf (R), 

± f f{x)^dx) = \f f dxf{x) " d y f{y \ t {dx)^{dy)- \ f V'(x)f'(x)n t (dx). 

al JR L J JR 2 x — V L JR 

This also proves the existence of weak solutions of the McKean-Vlasov equation. If one can 
further prove the uniqueness of the weak solution to the McKean-Vlasov equation, we can 
then derive that the empirical measure Ljy(t) converges weakly to /j,(t) as N — > oo. 

1.3 Interacting particle system 

The GDBM is a special example of interacting particle systems which can be defined by a 
stochastic differential equation of the form 

dx\ = dB\ - VV(xi)dt - ^J2^ W ^t - ^dt, i = l,...,N, 

where V is a external potential, W is an interacting function between particles, and B\ are 
i.i.d. Brownian motions on R d . See [37, 5, 11, 27, 39, 40] and reference therein. For N = oo, 
see e.g. [22, 31]. Let 

1 N 

i=l 

be the empirical measure of the particle system {xi(t),i = 1, . . . , N}. Then the above SDE 
can be rewritten as 

dx\ = dB\ - VV(xi)dt - W(W * n?){x\)dt. 

Under various assumptions which require that V and W are Lipschitz, as N — >• oo, it has 
been proved that ^ converges weakly to a measure- valued process fit on R, i.e., 
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which is the law of a nonlinear diffusion process on R defined by 

dX t = dB t - VV(X t )dt - V(W * (H)(X t )dt. 

where B t is a Brownian motion on R d . See [27]. Suppose that fj, t « dx, then the density 
function u = satisfies the nonlinear McKean-Vlasov equation (called also the nonlinear 
Fokkcr-Planck equation in [5] etc.) 

d t u = div(uV(logu + V + W *u)). (3) 

In [5], Benedetto, Caglioti, Carrillo and Pulvirenti developed the i 1 -theory of the non- 
linear McKean-Vlasov equation (3) in one-dimensional granular media. Assuming that the 
interaction function W and the potential V are Lipschitz and convex functions, they proved 
that the free energy functional 

F(u) = / u log udx + / Vudx + / / W (x — y)u(x)u(y)dxdy 
Jm Jm J Jm 2 

has a unique minimum u^, and it holds that 

\\ut — \\ L i — > as t — >• oo. 

Moreover, under the Lipschitz condition on V and W, they proved that there exists a 
constant C > such that 

W^uvt) < e c 'Wi(Mo,^o), 

where /j, t and v t are solutions to the nonlinear McKean-Vlasov equation (3) with initial dates 
and vq respectively. See also [37, 23] and references therein. 

However, the logarithmic Coulomb interaction function appeared in the distribution of 
the eigenvalues of /3- invariant random matrices ensemble, i.e., 

W{x) = - log|ar| 

is not a Lipschitz function on R. Thus, it does not satisfy the Lipschitz condition required 
in [27, 5], and the L 1 -theory of Benedetto et al.[5] does not apply directly to (3) associ- 
ated to the GDBM, even in the case of Gaussian ensembles, i.e., V(x) = As far as we 
know, the nonlinear McKean-Vlasov equation (3) with the logarithmic Coulomb interaction 
W(x) = — log \x\ and general external potential V has not been well studied in the litera- 
ture. In particular, even though the existence of weak solution of (3) can be easily derived 
from the McKean-Vlasov limit of any convergent subsequence of the empirical measure of 
the generalized Dyson Brownian motion, the problem of the uniqueness and the longtime 
asymptotic behavior of the solutions of the nonlinear McKean-Vlasov equation (3) remain 
as an open problem in the literature. 

1.4 Otto's calculus and the gradient flow of the Voiculescu entropy 

In [11], Carrillo, McCann and Villani used Otto's infinite dimensional differential calculus 
[32] on the Wasserstein space over R™ to study the convergence rate problem of the nonlinear 
McKean-Vlasov equation (3) in the granular media. Under some growth and smoothness 
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assumptions of the interaction potential W and the external potential V, they proved that, 
if V 2 W > and V 2 V > A, or if V 2 W > A and V 2 V > 0, where A > is a constant, then 
the solution u(t) of (3) converges to the equilibrium in the WVWassersten distance with 
an exponential rate, i.e., 

W 2 (u t , Uoo ) =0(e- xt ). 
Note that, for W{x) = — log |x|, we have 

V 2 W(x) = -^, x^O. 

Thus, the logarithmic Coulomb interaction potential has a strong convexity near its sin- 
gularity point x = 0. This suggests us to adopt the infinite dimensional calculus on the 
Wasserstein space initiated by Otto [32] and developed by Carrillo, McCann and Villani [11] 
to study the uniqueness problem and the long time asymptotic behavior of the solution of 
the McKean-Vlasov evolution equation. 

Let V be a C 1 function on R. According to Voiculescu [41], we introduce the free entropy 
as follows 

s vO) = - / / log \x - y\d/j,(x)dn(y) + / V(x)dfj,(x). 

J JR 2 JR 

It has the following electrostatic interpretation: Suppose that electrons are distributed on 
the one-dimensional space R in the presence of an external field with potential V. For a 
probability distribution /i of electrons, the electrostatic repulsion is given by 




log \x- y\ 1 dfi(x)dfj,(y) 7 



and the energy from the external field is J R V(x)dfi(x). Thus, the Voicuselecu free entropy 
is the total energy of electrons in an external field with potential V. When V satisfies 
the growth condition (1), it is shown in [20] that there exists a unique minimizer of the 
Voiculescu free entropy 

(i v = argmin AJeP(R) Ey(/i), 
called the equilibrium distribution of Ey. The relative free entropy is defined as follows 

Z V (n\(l V ) = Ey(/i) - Sy(/iy). 

Following Voiculescu [41] and Biane [3] , the relative free Fisher information is defined as 
follows 

where 

Jrx-v 
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is the Hilbcrt transform of /i. 

Note that, as the equilibrium measure jiy satisfies the equation 




we have 



Ivifiv) = 0. 

Inspired by Otto [32], Carrillo-McCann-Villani [11], Villani [39, 40] and Biane [3], we 
have the following result, which has been known by Biane and Speicher [4]. 

Theorem 1.1 (Biane and Speicher [4]) Let V <G C 2 (M). Then the McKean-Vlasov equation 
(5) is the gradient flow of the Voiculescu free entropy £y on the Wasserstein space ^(R) 
over R equipped with Otto 's infinite dimensional Riemannian structure. 

For the definition of the Otto's infinite dimensional Riemannian structure on the Wasser- 
stein space ^(R), see Section 4. 

1.5 Main results 

We now state the main results of this paper. Our first result gives the existence and unique- 
ness of the strong solution to the SDE of the generalized Dyson Brownian motion under 
reasonable condition on the external potential. 

Theorem 1.2 Let (W 1 , . . . , W N ) be an N -dimensional Brownian motion in a probability 
space (0,P) equipped with a filtration & — {JPt,t > 0}. Let V € C 2 (M) be a function 
satisfying the growth condition (1) and the following conditions 
(i) For all R > 0, there is Kji > 0, such that for all x, y € R with \x\, \y\ < R, 



Then, for any (3 > 1, and for any given Xn(0) = (Aj^(O), . . . , Ajy(0)) S An, there exists a 
unique strong solution (Ajv(£))t>o with infinite lifetime to SDE (2), with initial value Aat(0) 
and such that Ajv(i) € for all t > 0. 

The second result of this paper is the following result concerning the existence and 
uniqueness of the weak solution to the McKean-Vlasov equation, as well as the convergence 
of the empirical measure L^{t) towards the solution of the McKean-Vlasov equation. 

Theorem 1.3 (i) Suppose that V be a C 2 function satisfying the same condition as in The- 
orem 1.2. Suppose that 



(x-y)(V'(x)-V'(y))>-K R \x 




(ii) For some 7 > 0, such that 



-xV'(x) < 7(1 + |a;| 2 ), VxeR 



(4) 
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1 N 

Ln ^ = «E ^(0) -> A* e £»(R) as TV oo. 



TV 
fe=i 



Then, the family {LN(t),t G [0,T]]} is tight, and the limit of any weakly convergent subse- 
quence of {LN(t),t G [0, T]]} is a weak solution of the McKean-Vlasov equation, i.e., for all 

/eC^R), te [0,T], 

4 / f(x)Mdx)= 1 - I I MW^MVI {dx )*(dv)-± I V'(x)f'(x)p, t (dx). (5) 
dt J R 2J j R ! x-y 2 J R 

Equivalently, the probability density of fi t with respect to the Lebesgue measure on R satisfies 
the nonlinear Fokker-Planck equation 

d *> "(JIv-e^Y ( 6) 



dt dxV\2 

where 



Rp(x) = P.V. f ^^dy 
Jwx-y 



is the Hilbert transform of p. 

(ii) Suppose that there exists a constant K £ R such that 

V"(x) >K, lei. 

Let Hi(t) be two solutions of the McKean-Vlasov equation (5) with initial datas p,i(0), i = 1,2. 
Then for all t > 0, we have 

W 2 (Aii(t),M2(t)) < e- Kt W 2 (M0)^2(0)). 
In particular, the McKean-Vlasov equation (5) has a unique solution. 

(Hi) Let V be a C 2 function satisfying the same condition as in Theorem 1.2 and V" > K 
for some constant A"gl. Then the empirical measure L^it) weakly converges to the unique 
solution p{t) of the McKean-Vlasov equation (5). 

The following result gives us the longtime asymptotic behavior and the convergence rate 
of the weak solution of the McKean-Vlasov equation to the equilibrium measure of the 
Voiculescu free entropy for C 2 -convex potentials. 

Theorem 1.4 (i) Let V be a C 2 -convex potential. Then p{t) converges to py with respect 
to the Wasserstein distance in ^(R); i-e., 

W2(fi(t), nv) -> as t -> oo. 

(ii) Suppose that there exists a constant A"gl such that 

V"(x) >K, Vie R. 

Then for all t > 0, we have 

< e~ 2Kt £y(^ |Mv), 
W 2 {p u p v ) < e- Kt W 2 (p ,Hv)- 



9 



In particular, if K > 0, then p t converges to py w ith respect to the Wi-Wasserstein distance 
with the exponential rate K . 

(in) Suppose that V is a C 2 , convex and there exists a constant r > such that 

V"{x) >K>0, \x\> r. 

Then there exist constants C\ > and C2 > such that p(t) converges to \iy with respect 
to the Wi-Wasser stein distance with an exponential rate 

e -C!t 

W 2 2 (Mt,M < 5V(At |A*). t>0. 

1.6 Remarks 

Remark 1.5 (i) In [36], Rogers and Shi proved (even though did not state) the existence 
and uniqueness of the generalized Dyson Brownian motion under the condition that the 
potential V satisfies 

-xV'(x) < 7, Viet (7) 

It is clear that our condition (4) on V is weaker than Rogers-Shi's condition (7). To prove 
Theorem 1.2, we need to use the non-explosion criterion of the Bessell type SDE on the half 
line R + (cf. [19] p. 235-237), which was kindly pointed to us by Yves Le Jan in June 2011. 
(ii) The conclusion of Theorem 1.6 says that, for any given Ajy(O) < A^(0) < ... < 
A^ _1 (0) < \%(0), SDE (2) admits a unique strong solution with infinite lifetime such that 
X]^(t) < \ N (t) < ... < A^ _1 (t) < Aj^(t). Therefore, the generalized Dyson Brownian 
particle system (A]y(t), . . . , Aj^(t)) does not self-intersect for all time t > 0. 

Remark 1.6 Let [3 = 1,2,4. Let V be an analytic function satisfying the assumptions of 
Theorem 1.2. Following an original idea due to Dyson [15] and developed by other authors 
[12, 36, 2, 18] for the special case V{x) = 9x 2 , we can prove that the generalized Dyson 
Brownian motion Xj^(t) < . . . < \%{t) can be realized as the eigenvalues of a matrix-valued 
diffusion process X t , which satisfies the following SDE 

dXt = \[j^ dW t - \nrV{X t )dt, (8) 

where B t is the standard Brownian motion on VP N (see Section 2.1 ), and V denotes the gradi- 
ent operator on %%. Moreover, Ajv(i) = (Ai(t), . . . , Ajv(i)) * s a An -valued semi-martingale. 
See Theorem 2.2 and Theorem 2.3. 

Remark 1.7 In [3], Biane pointed out that there exists non- convex potentials such that the 
equation Hpi(x) = ^V'(x) on the supports of fii holds for distinct measures fii — pi(x)dx. 
The simplest example due to Biane [3] is a two-well potential V satisfying V(x) = \{x — ai) 2 
in [di — 2,ai + 2] for points ai,a 2 with \a\ — a 2 | > 4, then the semi-circular measures centered 
on a\ and ai respectively satisfy Hpi(x) = V'(x) on their supports respectively. By Theorem 
1.2 (due to Biane and Speicher [4]), the McKean-Vlasov limit n{t) — pt{x)dx of Ljq(t) 
satisfies the free Fokker-Planck equation 
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It is natural to ask the question whether fi(t) converges to fly in the weak convergence 
topology or with respect to the W2-Wasserstein distance (or the W p -Wasserstein distance 
for some p > 1) for non-convex potentials V. In [4], it was pointed out that this cannot 
be true in general. Indeed, as /it satisfies the gradient flow of the Voiculescu free entropy 
Ey on ^(R), /j, t may converge to a local minimizer ofSy which is not necessary the global 
minimizer fiy. See also Section 6 and in particular Conjecture 6.1. 

The rest of this paper is organized as follows. In Section 2, we prove Theorem 1.2. In 
Section 3, we prove Theorem 1.3 (i). In Section 4, we prove Theorem 1.4. In Section 5, we 
prove Theorem 1.3 (ii) and (iii). Finally, we study the double- well potentials and raise some 
conjectures in Section 6. 

Acknowledgement. During the preparation of this paper, the authors have benefited 
from some discussions with G. Akemann, F. Gotze, K. Kuwae, Y. Le Jan, Y. Liu, T. Lyons, 
S.-G. Peng, M. Rockncr, W. Sun and M. Vcnker. We would like to express our gratitude to 
all of them for their interests and for useful comments. 

2 Generalized Dyson's Brownian motion (GDBM) 
2.1 Proof of Theorem 1.2 

For T e R+, we denote by C([0,T],R N ) the space of continuous process from [0,T] to R N 
and 0>(C([O,T],A N )) is the set of all probability measures on C([0, T], A N ). 

(1) Fix R > 0. Define the truncated Dyson Brownian motion by 

= J^dWl + 1 £ MW) \UWt - lv'(XUt))dt, (9) 

with \if tR (0) = A 4 Ar (0) for 1 < i < N, where <f> R (x) = x" 1 if \x\ > R' 1 , and </> R (x) = R 2 x 
if \x\ < R~ x . By Theorem 3.1.1 in [35], since <j)R is uniformly Lipschitz and V satisfies the 
conditions (i) and (ii), SDE (9) has a unique strong solution, adapted to the filtration & ' . 
Let 

t r := inf{t : min | \^ R (t) - X NM (t) \< R- 1 }. 

Then t r is monotone increasing in R and 

X N>R (t) = \ N ,R>(t) for all t<r R and R < R' . (10) 

(2) We now construct a solution to SDE (2) by taking Ajv(t) = Ajv,ii(i) on the event 
[tr >t] = {\\%(s) - A J w (s)| > R- 1 , Vs < t,\/i ^ j}. To obtain a global solution {X N (t),t € 
M+} to SDE (2), we need only to prove that t r tends to infinity almost surely when R — > oo. 

To this end, let us consider the Lyapounov function 

1 N 1 

f(xi,...,X N ) = Ji^V(Xi) - j^^\0g\Xi-Xj\. (11) 

1=1 ijij 



11 



Using log \x — y\ < log(|a;| + 1) + log(|y| + 1) and (1), we can derive that 

/(a:i,...,a;j V )>-(l + J)Iog2. (12) 
and V(x) - 21og(|x| + 1) > C for C := -2(1 + S) log 2. Moreover, for all i ^ j, it holds 



jplog\Xi-Xj\ < f(X!,...,X N )-C. (13) 



For any M > 0, set 

T M = inf{t > : f(X N (t)) > M}. 
Then, Tm is a stopping time. Moreover, on {Tm > T}, we have 

lAjv-(t) - A^(t)| > JT 1 := e N2 (- M+c \ V i < T. 

Thus, on the event [T < Tm], {^N,ii(t), t < T} is an adapted solution to SDE (2). It remains 
to prove that for all t > 0, we have 

¥(3M e N : T M >t) = l. (14) 

(3) To prove (14), we need only to prove that, for all K > 0, we have 

P(3M e N : T M A Or > t) = 1, (15) 

where 

Or = inf{i > : Aiv(i) ^ [-if, X], for some i = 1, . . . , N}. 

By (11) and (13), to prove (15), it is equivalent to show that almost surely Ajy(i) and X J N (t) 
never collide up to (m- We shall prove this claim below. 

(4) By Ito's formula, for all / e C 2 (K), we have 



1 N ( 1 1 

+2^E(-i i "(a-„w)I 2 + ^E 

i=i y j^i 



A^W-A^W 



+7777^ I V'fAVft)) + 4 Y — — I dt + dM N (t), 

where Mn is the following local martingale 

dM N (t) = Y ( V'(\Ut)) - — V - , r 1 rT7 I <*W7 ■ 

By [18], we have 



N 

E 

i=l 



= 0. 
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Thus 

df(X N (t)) 



TV 



(it 



= dM N (t) + A N (t)dt, 
where, as /3 > 1, it holds 



1 / 1 N 3 



y'(Aj v (t))-r(A^(t)) 



AVW-A^(t) 

We now prove the following lemma, which is a stronger version of the claim in (3). 

Lemma 2.1 Suppose that the processes (Ajv(t))t>o defined by (2), at least up to the stopping 
time 

C = inf{i : Aiv(t) - A^(i) 3^jor AV(t) - oo 3 i}. 

Then 

P(C = oo) = 1. 

Proof. Fix X > 0, T > and P > such that A] v (0) e [-if, if] for all i = 1,...,N, 
and |AV(0) - Aj vr (0)| > P" 1 for all i ^ j, i,j = 1,...,N. Let Ci(if) > be such that 
V"(x) < d(K) for all x e [-if, A]. Then, < Ci(if) and {/(Ajv(i A Or ) - C X {K )(t A 

Cat), f € [0,T]} is a super-martingale. 

Let A := {r fl < Cx,^ < T}, and C 2 (if) := inf{V(x) : |x| < if}. Then 

/(Ajv(0)) + TC^K) > E(f(X N (T A( K A r R ))) 
= E (f(\ N (T R ))l A ) + E (/(Ajv(T A Ck))Uc) 

> (JL log R - -L (TV 2 - TV - 1) log(2if ) + C 2 (J0) P(A) 
-^(7V 2 - AT) log(2P) + C 2 (K)j F(A C ) 
± logP + log(2if)) P(A) - ^_(7V 2 - AT) log(2tf) + C 2 (tf ), 
whence 

p, . , ^ TW ^(/(Aat (0)) + TCi(ff)) + N(N - 1) log(2g) - C 2 (iQ 
P(r* < CK,r fl < T) < _____ . 

Taking R -> oo, for all if and T, we have P(too < £k A T) = 0. Letting T tend to infinity, 
we obtain P(too < (k) — 0, which yields that ^ log \X z N {t) — X J N (t)\ > — oo almost surely 

for all t < Qk- Moreover, letting K tend to infinity, we get 

F(Too < C) - 0, 
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where ( := lim (k = inf{t > : /(Ajv(i)) = 00} is the explosion time of /(Ajv(i))- This 

K— >oo 

means that the particles Ajy-(t), . . . , A^(i), could not collide before the explosion. 
We now prove that £ is infinity almost surely. To this end, let 

1 N 

where /(a;) = Now ■ A"^*- = By Ito's formula, we have 



dRt 



j=l 

By introducing a new Brownian motion £?, we have 



Under the assumption (4), and using the comparison theorem of one dimensional SDEs, cf. 
[19], we can derive that 

R t <R' t , Vi>0, a.s, 

where R' solves 

2 jyf / ^ TV 1 1 \ 

dR * = N^i m + {f3N + ^r + r< + ^) *■ 

with R' = 0. Moreover, by Example 8.2, (8.12) in Ikeda and Watanabe [19] (p. 235-237), 
the process R' never explodes. So the process R cannot explode in finite time. 

(5) By the continuity of the trjectory of Ajv(i)> it is easy to see that Ajv(t) G Ajy for all 
t > 0. The proof of Theorem 1.2 is completed. □ 



2.2 From matrix diffusions to GDBM 

The following result indicates a way to introduce generalized Dyson Brownian motion as the 
eigenvalues process of a matrix valued diffusion process. To simplify notation, let VP N be 
the ensemble of N x N matrices: for f3 = 1, it denotes the N x N real symmetric ensemble, 
for (3 = 2, it denotes the N x N Hermitian complex ensemble, and for (3 = 4, it denotes the 
N x N symplectic ensemble. 

Theorem 2.2 Let (3 = 1,2,4, and V : R — »■ M &e a reaZ analytic function. Let X t be a 
VP N -valued diffusion process defined by 

dX t = J^dB t - l -VTiV{X t )dt, (16) 
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where B t is the standard Brownian motion on T-L^ N , V denotes the gradient operator on W^. 
Let Ajv(t) = (^jv(i)) ■ • • > Ajv(O) be eigenvalues ofX^it)- Then, Ajv(t) satisfies the SDE 
(2) /or i/ie generalized Dyson Brownian motion with (3 = 1,2, 4. 

Proof. We only prove Theorem 2.2 for = 2. Note that, for analytic function V(x) — 

oo 

a kX k on R, we have 

fe=0 

oo 

VWpf) = J2 ka k X k - x = V'(X). 
fe=i 

Hence the SDE (16) for X t can be written as follows 

dX t = ^=dBt - \V(X)dt. (17) 

Let X\{t) < . . . < XN{t) be the ordered eigenvalues of X t . Let 

/ : X -> D = diag(xi, . . .,x N ) 

be the matrix transformation such that X = UDU* , where U — (u\, . . . ,Un) is an unitary 
matrix. Equivalently, we have 

i.e., Ui is the eigenvector of M with eigenvalue x t . Write / = (/i, . . . , fjy), where fi(X) = x i} 
i = 1, . . . , N. By Ito's formula, we have 

dxi(t) = V dx Ji(X t ) + l -V 2 fi{X t ){dX u dX t ). (18) 

By the first order Hadamard variational formula, see [38] , we have 
VdxJi(X) = u*dX tUi 

= -±=u*dB tUi - ^u*V'(X t ) Ui dt. 

Note that 

u*V'(X t ) Ui = V'( Xi (t)). 

By the rotational invariance of the Brownian motion, U*B t U is a Brownian motion on C N . 
Denote 

W l t =u*B t Ui, i = l,...,N. 

Then 

V dx Ji{X) = -^=dWi - l -V\x t {t))dt. (19) 



On the other hand, by the second order Hadamard variational formula, see [38], it holds 

\u*dX t Ui\ 2 
W 2 h{X t ){dX t ,dX t ) = 2J2 3 



_ y ii(i)-a!j(t)' 
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By Ito's calculus, we have 



1 *WR 1 
— UjtilhlU--,^ 

1 . ,„ 1 



\u*dX tUl \ 2 = \—u*dB tUi - -u*V'(X t ) Ul \ 2 



= l-^u^dBtUt-^V'ix^t^udtl 2 



^=u*dB t u^ 2 
N 



where we have used the fact that U*B t U is a Brownian motion on C N . Hence 

^MX t )idX u dX t ) = ^ ^ W) dt. (20) 

From (18), (19) and (20), we derive that (xi(t), . . . , xjv(O) satisfies the following SDE 

= ±™n - \vww + ^ E ^y^)^- ( 21 ) 

The proof of Theorem 2.2 is completed. □ 
2.3 Prom GDBM to matrix diffusions 

The following result provides a random matrix representation for the generalized Dyson 
Brownian motion which is defined by solving SDE (2). 

Theorem 2.3 Let (3 = 1,2 and Xn{0) G An- Then, there exists a Sn (respectively, Wn)- 
valued diffusion process (X N 'P(t),t > 0), such that its eigenvalues process (Ajv(i))t>o is a 
solution of the SDE (2) for the generalized Dyson Brownian motion. 

Proof. We only prove Theorem 2.3 for (3 — 1. Without loss of generality, we assume 
that X N (0) is the diagonal matrix D=diag (Ajy, . . . , Ajy). Let M > be fixed. We consider 
the strong solution Ajv(i) of SDE (2) untill the stopping time Tm- We let Wij, 1 < i < j < N 
be independent Brownian motions. Hereafter, all solutions will be equipped with the natural 
filtration Tt — cr(wij(s),Wi(s), s < t A Tm), where Wi the Brownian motions of SDE (2), 
independent of Wij, 1 < i < j < N. For i < j, define Rij(t) by solving SDE 

dRfJt) = — dwijtt), R$(0) = 0. (22) 

For i > j, set R^(t) = -Rf t (t). Let N be the strong solution of SDE 1 

\o N {t)d{{R N f,R N ) u 



dO N (t) = N (t)dR N (t) - l-0 N (t)d((R N ) T ,R N ) t , 



O N (0) = I 



N 

1 Here, for two semi-martingales A and B with values in Sjv, (A, B)t = ( J2 (A%k> ^fcj)t)l<» j<iv is the 

fc=i 

martingale bracket of A and B, and (A)t is the finite variation part of A at time t. 
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Since SDE (23) has uniformly Lipschitz coefficients, we obtain the existence and uniqueness 
of strong solutions of (23) with respect to the filtration T t up to the stopping time Tm- By 
Lemma 4.3.4 in [2], we have N (t) T N {t) = I at all times. 
Let Y N (t) = N (t) T D(X N (t))0 N (t) and define W N (t) by 

dW N (t) = (0 N (t)) T dY N (t)0 N (t) 

with the initial values W N (0) = Y N (0) = X N (0). 

By the same argument as used in [2, 18], we can prove that 

dWi$(t) = \f^dW l t - ±V'(X N (t))dt, V* = 1, . . . , N, (24) 

and for i ^ j, 

dWiZ(t) = ] p^dw ij (t). (25) 

Since (0 N (t),t > 0) is adapted, dY N (t) = N (t)dW N (t)(0 N (t)) T is a continuous matrix- 
valued semi-martingale. Set 

~ jwij(t), iii^j, 
W t = { 

\\V l (t), Hi=j. 

Then we can write dY N (t) = O n (t)dW N (t)(0 N (t)) T as follows 



dY N (t) = N (t)y^dW t -±D(V\\ N (t))dtj(0 N (t)) T 

= \[^0 N (t)dW t {0 N {t)) T l -0 N (t)D{V'(\ N (t)){0 N (t)) T dt. 



Set 



Bt= I N ( S )dW s (0 N (s)f. 



By Levy's characterization of Brownian motion, B t is a Sjv-valued Brownian motion with 
respect to Ft = cr(wij(s), Wi(s), s < t). Moreover, by the fact V is a real analytical function 
and Y N (t) = N {t)D(\ N (t))(0 N (t)) T 1 we have 

VTvy(y t ) - o N (t)D(V(\ N (t))(o N (t)) T . 

Hence 

dY N {t) = \f^dB t - ^VTrV (Y t )dt. (26) 

Thus, ljv(i) is the Sjv-valued diffusion process with generator L — ^^s N — iVVTrV-V. Note 
that the SDE (26) is as the same type as SDE (16) that X N (t) satisfies, and Y N (0) = X N (0). 
By the uniqueness of weak solution to SDE, ljv(t) has the same law as Xj^(t), and Ajv(i) is 
the eigenvalue process of Yjv(£)- The proof of Theorem 2.3 is completed. □ 
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Remark 2.4 Theorem 2.3 can be considered as a reverse of Theorem 2.2. The way from 
Xjv(t) to Ajv(t) in Theorem 2.2 is more direct and does not need to deal with the question 
of the existence and uniqueness of strong solution of the SDE (2) with singular drift (i.e., 
Theorem 1.2). On the other hand, the reverse way from Ajv(t) to Y)v(t) in Theorem 2.3 
shows that the generalized Dyson Brownian motion Ajv(t) must be obtained by the way from 
Xpf(t) to Ajv(t) as we did in Theorem 2.2. Moreover, Theorem 2.3 also shows that the SDE 
(2) for generalized Dyson Brownian motion has the uniqueness in distribution. We refer the 
reader to [30, 18, 2, 8] and references therein for further work on matrix-valued diffusion 
processes. 



2.4 Ito's calculus 

Let (W 1 , . . . , W N ) be independent Brownian motions and (Aj v (0), . . . , Ajy(O)) be real num- 
bers, let f3 > 1 and let \N(t) t >o be the unique strong solution to SDE (2). Then by Ito's 
calculus, we know that for all / € C 2 ([0,T] x M,R), 



f f{t,x)L N (t,dx) = [ f(0,x)L N (0,dx)+ [ f d s f( 

Jr Jr JO JS, 



s, x)Ln(s, dx)ds 



Iff / d x f(s,x)-d y f(s,y) 



o J Jr 2 x ~y 
t 



-L N (s, dx)L N (s, dy)ds 



T7 / / dlf(s,x)L N (s,dx)ds 
JV Jo Jr 

I f fv'(x)d x f(s,x)L N (s,dx)ds + M t N (t), 



JR 



where M* N is the martingale given by 



1 / 2 N rt 



i=i J ° 



M N^ = if\li^y. I d x f{s,y N { s ))dwi, vt<r. 



N\l f3N 
Note that 



(3N 2 



3 McKean-Vlasov limit as N — > oo 
3.1 Proof of Theorem 1.3 (i) 

We first prove the tightness of {L N (t),t e [0,T]} in C([0,T], &>(R)), then we show that 
the limit of any convergent subsequence of {Ljv(t),£ € [O,? 1 ]} satisfies the McKean-Vlasov 
equation (5). Here, for T e M + , we denote by C([0, T], ^(R)) the space of continuous 
processes from [0,T] into ^(R) equipped with the weak convergence topology. 

3.1.1 The tightness of the empirical measures 
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We follow the argument used in [36] to prove the tightness of {Ljv(f),i € [0,T]}. Let us 
pick functions fj G C b °°(IR, C),j = 1,2, ... , which is dense in C 6 (R). Thus 

= (m'./j), Vj n = n'. 

We also pick a C°° function / : R — > [1, oo) with the properties 

fo(x) = fo(-x), f (x) ^ oo as noo, 

Taking test functions in the Schwartz class of smooth functions whose derivatives (up to 
second order) are rapidly decreasing, we may assume that 

fj, f", V'f'j arc uniformly bounded for all j > 1. 

By Ethier and Kurtz [16] (p.107), to prove the tightness of {L N (t), t G [0,T], N > 1}, it is 
sufficient to prove that for each j the sequence of continuous real-valued functions 

{(L N (t),fj), te [0,T], iV> 1} 

is tight. To this end, note that, by Theorem 1.2, there is non-collision and non-explosion for 
the particles X l N (t) for all t G [0, oo). By Ito's formula in Section 2.2, we have 

d{L N (t)J) = l^^/'^t)^ 

fl{ - x) - fl(y) L N (t,dx)L N {t,dy)dt. 
2 J Jr2 a; - y 

This yields 

(L N (t),fj) = (MO),/,)!- / / jV ^ ^(a.^Ljy^dy)^ 

z Jo J JR 2 x — y 

~\ J\l n {8),V 'fyds + J\l n {s), (j - lj ^fj)ds + M% (t) 
= h(N)+I 2 (N) + I 3 (N)+I 4 (N) + M%(t), (27) 

where 

i ht - /•* N 

Note that, as £jv(0) is weakly convergent, h(N) is convergent. By the assumption that fj 
and fj are uniformly bounded (hence /j are uniformly bounded) , we can easily show that 

{M^{t),t G [0, T]} and Ii{N) converge to zero. Moreover, by the assumption that V'f'j 
and fj are uniformly bounded, the Arzela-Ascoli theorem implies that I2(N) and Is(N) 
are tight in C([0,T],R). Thus the sequence {(Ljv(*))t>o : N > 1} is tight in C([0,T],R). 
Tightness also follows for j = if we have 

(.Ljv(O), /o) — > finite limit as TV — > oo. 
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So let us suppose that the initial distribution L^(0) have the property (Ljv(O), fa) < K for 
some K, for all N. For given we could always find Ljv(O) and /o to satisfy this and the 
other conditions, and this gives the tightness for j = also. 

3.1.2 Identifying of the limit process 

Without loss of generality, assuming that {LN j (t),t e [0, T}} is a convergent subsequence 
in C([0,T], ^(R)). Then, for all / € Cf(R), the Ito's formula (27) in Section 2.2 and the 
argument used in Section 3.1 show that = lim (Ln f) satisfies the following 

j->oo 

equation 

f f(x)fH(dx) = [ f(x)H>{dx) + I f ( [ dxf{x) ~ dyf{y) tis(dx)v s (dy)ds 
Jr Jr 2 J J J R2 x-y 

~\ I [ V'{x)f'{x) l x s {dx)ds. 
L Jo Jr 

That is to say, n{t) is a weak solution to the McKean-Vlasov equation (5). Suppose that \x t 
is absolutely continuous with respect to the Lebesgue measure dx, and denote 

Pt{x) = ^ 

Integrating by parts in (5), we can prove that p t satisfies the nonlinear Fokker-Planck equa- 
tion (6). The proof of Theorem 1.3 (i) is completed. □ 

Remark 3.1 To characterize the McKean-Vlasov limit \i t of the family of empirical mea- 
sures {Lj^(t),t G [0,oo)} 7 we need only to use the test function f(x) = (z — x)^ 1 , where 
z E C\M, instead of using all test functions f € C 2 (R) in the McKean- Vlasov equation (5) . 
Let 

GAz) = > 



«(*) = / « 
Jr z 



x 

be the Cauchy-Stieltjes-Hilbert transform of \x t . Then G t (z) satisfies the following equation 

d d 1 ( V'(x) 

-G t (z) = -Gt{z) — G t {z) - - ^ j-U^^dx). 

In particular, in the case V{x) = 8x 2 , since 



J 

Jr 



x d 

Ht{dx) = z—G t (z) + G t {z), 



(z — x) 2 dz 
we obtain 

-G t (z) = (-G t (z) + dz) Q-G t {z) + 8G t (z). 

Remark 3.2 Following [18, 2], we can also prove a stronger version of the tightness of the 
empirical measure. More precisely, let T e R+, and assuming V satisfy the conditions in 
Theorem 1.6. Assume that 



sup / log(x + 1)Ln(0, dx) < +oo. 

JVeN. 
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Then, for all T e R+ and L e N, there exists a compact set K(L) C C([0,T], ^(R)) sucft 
that 

¥(L N (-) € K(L) C ) < e- N " L . 

In particular, the law of (Ljv(s),s € [0, T]) is almost surely tight in C([0,T], <^(R)). XTiis 
w/Z 6e used m our forthcoming paper for a study of the large deviation principle for the 
generalized Dyson Brownian motion. For a proof see our forthcoming paper. 



4 McKean-Vlasov equation and optimal transport 

In Section 3.1, we have proved Theorem 1.3 (i), which asserts the existence of weak solution 
to the McKean-Vlasov equation (5), equivalently, the existence of weak solution of the non- 
linear Fokker-Planck equation (6). In this section, we use the mass optimal transportation 
theory to study the uniqueness and the longtime asymptotic behavior of the McKean-Vlasov 
equation (5) and the nonlinear Fokker-Planck equation (6). 
Let 

W(x) = — 21og|ar|, x ^ 0. 

Then the nonlinear Fokker-Planck equation (6) can be rewritten as follows 

d tP = V -( P V(V + W *p)). (28) 

Before going to study the nonlinear Fokker-Planck equation (6) (i.e., (28)), we first recall 
some results due to Carrillo, McCann and Villani. In [11], Carrillo, McCann and Villani 
studied the McKean-Vlasov evolution equation of the granular media, which is given by the 
following 

dtp = V- (pV(logp + V + W*p)). (29) 

They proved that the McKean-Vlasov equation can be realized as a gradient flow of a free 
energy functional on the infinite Wasserstein space. More precisely, we have 

Theorem 4.1 (Carillo-McCann-Villani[ll]) Let 

F(f)= [ plogpdv+ [ pVdv + \[ f W(x-y)p(x)p(y)dxdy. (30) 

Then the McKean-Vlasov equation (29) is the gradient flow of F with respect to the following 
infinite dimensional Riemannian metric on ^(R^) (cf. Otto [32]): 

9fdv{si,s 2 ) = / sxs 2 fdv, 

JR d 

where fdv G V 2 (R d ), Si, s 2 € T fdv ^> 2 (^ d ) = {s : R d -> M : J M sdv = 0}, and 

s, = -V.(/Vj>j) 

for some Pl G W^ 2 {R d ), i = 1,2. 



21 



Moreover, based on Otto's infinite dimensional geometric calculation on the Wasserstein 
space, Carrillo, McCann and Villani [11] proved the following entropy dissipation formula 

Theorem 4.2 (Carrillo-McCann-Villani[ll]) Denote £ := V(logp + V + W * p). Then 

- / Pltfdv, (31) 

JR™ 

2 f pTr(D^ T (D^dx + 2 [ (D 2 V • pdx 

jr™ jr™ 

+ [ (D 2 W(x-y)-[Z(x)-Z(y)],[li(x)-li(y)})dp(x)dp(y). (32) 

JR2™ 

Now let V : K — > [0, oo) be a C 2 function with growth condition (1). In [4], Biane and 
Speicher proved that the free Fokker-Planck equation 

£~£<*<H*-in> 

is the gradient flow of £y on the Wasserstein space ^ 2 (R)' See also [8]. 

By analogue of the proof of Theorem 4.2 in [11] , and observing that for W(x) = —2 log \x\ 
it holds 

f= V(F + W * p) = V - 2Hp, 
we can derive the following dissipation formula for the Voiculescu free entropy. 
Theorem 4.3 Let £ = V - 2Hp. We have 

-2 f [V'{x) - 2{Hp){x)f p(x)dx, (33) 
Jr 

2 / v"{x)\V'{x) - 2Hp(x)\ 2 p(x)dx 
Jr 

+ / / J^m^ss^i^ ^^ ( 34, 

J Jr2 (a; - yy 

Proof of Theorem 1.4 (i). By Corollary 3.2 in Biane [3], for any convex V, there exists 
a unique equilibrium measure p (indeed p = pv) with a density p satisfying the Eulcr- 
Lagrange equation Hp(x) — ^V'(x) for all x <E supp(^i). Thus, £y has a unique minimizer 
pv- By the fact that £y is lower semi continuous with respect to the weak convergence 
topology, see e.g. [2, 18], we see that it is also lower semi continuous with respect to the 
Wasserstein topology on ^(R). Moreover, as V is C 2 -convex, Theorem 4.3 implies that £y 
is a (displacement) convex functional on & 2 (R) • 

By Proposition 4.1 in Klockncr [24], we know that .^(R) has vanishing sectional cur- 
vature in the sense of Alexandrov. More precisely, for any pi 7 p 2 ,p3 € £P 2 (R) and for any 
Wasserstein geodesic 7 : [0, 1] — >• ^ 2 (R) such that 7(0) = pi and 7(1) = p2, for all t € [0, 1], 
it holds that 

W%(p 3 ,l(t)) = tW 2 {p^pi) + (1 - t)Wi(p 3 ,p 2 ) - t(l - t)W 2 (p ll p 2 ). 



dt F ^ = 



dt 2 



npt) = 



—T, v (p t \pv) 



dt 2 



^viMpv) 
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Therefore, .^(R) is a nonpositively curved (NPC) space in the sense of Alcxandrov (even 
though ^ 2 (]R") is an Alexander space with nonnegative curvature for n > 2, see e.g. [1]). 
By Mayer [29], we can conclude that W2(p,t, P-v) — > holds if we only assume that V is a 
C 2 -convex potential. The proof of Theorem 1.4 (i) is completed. 

□ 

Proof of Theorem 1.4 (ii). The proof follows the same argument as used in [32, 33, 11]. 
We use the gradient flow of the Voiculescu entropy on the Wasserstein space and the 
A-convexity along the geodesic displacement between two probability measures. 
Since V" > K, we have 

d 2 

■5V(m*|mv) > K. 



alt 2 ' 



By Otto's calculus, we know that 
RessY, v (p t \p v ) 

which implies that 



0t ' dt J ~ dt 2 
HessSy(^) > A. 



Let p(0) = p(0)dx and p,(l) = p(l)dx be two probability measures with compact support 
on M, let p,(s) = p(s)dx be the unique geodesic in the Wasserstein space ^(R) linking /x(0) 
and p(l). Then 



d 2 
ds 



■j^Hvi^s)) = HessSy(p(s)) 



dp(s) dp(sY 
ds ' ds 



)>-* 


dp(s) 




ds 



Therefore, for some a* G (0, 1), 
S V (P(1)) - S v (p(0)) = 



> 



£Up(0)) + -££(p((r))| e 
/ rfp(o") 

^(g) 



(7 = 



if 



/* 


a P (<7) 


/o 


5a 



<j=0 



+ ^W 2 (p(0),p(l)). 



Similarly, 



S v (p(0))-Ev(p(l)) > 



gg(g) 



A 



W 2 2 (p(0),p(l)). 



<T = 1 



Summing the two inequalities together, we obtain 



/ Ma) 
\ da 



(7=1 



Ma) 
da 



> KW 2 (p(0),p(l)). 



<T = 



(35) 
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we 



Let pt(s,x)dx : [0, 1] — > ^2(K) be the unique geodesic between p, t and py. By Otto [32], 
have the following derivative formula of the Wasserstein distance 



dpt(s) 

ds 
dpt(s) 



dp t + 2 



3 = 



/ dpt(s) 
\ ds 



dpv 



ds 



(a;),VEv 



s=0 



(36) 



< -2KWi(p t7 p v ), 
where in the last step we have used (35). The Gronwall inequality yields 

Wi{p u pv)<e- 2Kt Wi{^^v). 

Recall that 

VSy(Mv) = 0. 

By the fact that p t is the gradient flow of Ey on ^(R) and using the uniform if-convexity 
of Ey , we have 

|||VMM*)llk(«) = 2^V||VEv(Mt)|lk(R),^] 
= -2HessJ: v (p t )l—,— 



-2K 



dfit 



dt 



Since Ey(^iy) = 0, we derive that 



-2K\\VX v (p t )\\% 2 



—Y> v {pt\pv) 



-HVEy^OIlk 
'°° d 



which implies 



f°° d 

= J -^\\VZ v {p s )\\% 2(R) ds 

/OO 
||VEv(/ia)|||. 2( R)ds 

f°° d 

= 2K —Z v (p s )ds 
= -2KX v (p t \p v ), 

Ey(/i t |Mv) < e- 2ift Ey(/i |/Uy)- 



□ 



The proof of Theorem 1.4 (ii) is completed 

To prove Theorem 1.4 (hi), we need the following free logarithmic Sobolev inequality 
and free Talagrand transportation cost inequality due to Ledoux and Popescu [25]. 
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Theorem 4.4 (Ledoux-Popescu [25]) Suppose that V is a C 2 , convex and there exists a 
constant r > such that 

V"(x) >K>0, \x\> r. 

Then there exists a constant c — C(K, r) > such that the free Log-Sobolev inequality holds: 
for all probability measure fi with Iv{p) < oo, 

2 

^vivW) < -Mm)- 

c 

Moreover, the free Talagrand transportation inequality holds: there exists a constant C = 
C(K, r, V) > such that 

CWifanv) <5V0«IM- 

Proof of Theorem 1.4 (Hi). By Biane and Speicher [4], we have the following entropy 
dissipation formula 

d 1 

— Y, v (n t \ij, v ) = --lv(jUt). 

By Theorem 4.4, there exists a constant C\ > such that the free LSI holds 

2 



which yields 



d C 

— E v (/z t |/ivO < — ^ s v(m*)- 



By the Gronwall inequality, we have 

Ey(^IM < e- Cl */ 4 S y (/x |^). 

By Theorem 4.4 again, there exists a constant C 2 > such that the free transportation cost 
inequality holds 

W£(fJ,,fi V ) < -^-E V (fJ,\fJ, V ). 

Therefore 

e -Cit/4 

W£{/J,t,Hv) < — ^ Ev(/i |A»v)- 

This finishes the proof of Theorem 1.3 (hi). □ 

5 Proof of Theorem 1.3 (ii) and (iii) 

Proof of Theorem 1.3 (ii). In the proof of Theorem 1.4, we have proved the following 
inequalities 



E, l)i2U)) - E, (/,.,.(/)) ( grad^Ey(p t (s)), -|-p t (s) 



+ ^Wi{p 1 {t),ii 2 {t)), 

s=0 Z 
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and 



M/*i(*)) - M/fc(*)) > - ( grad w S y (p t (s)), — p t (s) 



+ f W 2 2 ( Ml (i), M2 (t)). 



Summing them together, we obtain 



d 



gTad w Y, v (pt(s)), g^Pt(s) 



s=0 



grad w E y (p t (s)), ^Pt(s) 



<-W 2 2 ( M i(i),p 2 (i)). 



8=1 



By Otto [32], we have the following derivative formula of the Wasserstein distance 



-wi(Mt)^2(t)) 



dpt(s) 
ds 



(xUi 



dm(t) + 2 



s=0 



dpt(s) 



= 2 



ds 
dptjs) 
ds 



dp 2 (t) 



(x),VE v (/i2(t)) 



which implies 



W 2 (/xi(t),M*)) < e- A W 2 (/ii(0),/i2(0)) 



As a consequence, the McKean-Vlasov equation (5) has a unique weak solution. This finishes 
the proof of Theorem 1.3 (ii). □ 

Proof of Theorem 1.3 (Hi). By Theorem 1.3 (i), the family {L N (t),t <= [0,T]} is tight 
with respect to the weak convergence topology on &(M), and the limit of any convergent 
subsequence of {L^{t),t £ [0,T]} is a weak solution of (5). By the uniqueness of the weak 
solution to (5), we conclude that Ljv(i) weakly converges to p(t). The proof of Theorem 1.3 
(hi) is completed. □ 

By the same argument as used in Otto [32] and Otto-Villani [33], we can prove the 
following HWI inequalities. To save the length of the paper, we omit the proof. 

Theorem 5.1 (HWI inequalities) Suppose that there exists a constant K G R such that 

V"(x) >K, Vie R. 

Let pi e ^2(R), i = 1, 2. Then for all t > 0, we have 

K 



£ v (/ii) - E v (/z 2 ) < W A 2 (Mi,Ai2)||grad lv Ey(^ 1 )||^ 2(R) - —Wi{p 1 ,p 2 )- 
In particular, for any solution to the McKean-Vlasov equation (5), we have 

MM*)) < W 2 (Mt),M||grad w E y (Mt))||^ 2(K) ~ yW 2 2 (M*).M*))- 



(37) 



(38) 



where 



\gvad w T, v (p)\\% 2 



p\V'{x) - 2Hp(x)\ 2 dx. 
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6 Double-well potentials and some conjectures 



Theorem 1.3 and and Theorem 1.4 apply to V(x) = a\x\ p with a > and p > 2. When 
a = |, p = 2 and /3 = 1, 2, 4, it corresponds to the GUE, GOE and GSE. Moreover, Theorem 
1.3 and Theorem 1.4 also apply to the Kontsevich-Penner model on the Hcrmitian random 
matrices ensemble with external potential (cf. [10]) 

V(x) = ^ - ^- - c\og\x\. 

Note that, for all x ^ 0, 

V"(x) = ax 2 + ^-b 
x z 

> 2^-b>0 

provided that a > 0, c > and 4ac > b 2 . 
Let us consider the double well potential 

1 c 

4 2 

where c £ tt is a constant. By [20, 6], it has been known that the density function of the 
equilibrium measure fly is explicitly given as follows: 

(i) When c < -2, 

p(x) = ^-\x\\/{x 2 - a 2 )(b 2 - x 2 ) a < \x\ < 6, 
= otherwise, 

where a 2 = — 2 — c and b 2 = 2 — c. 

(m) When c = —2, = ^x 2 ^ — x 2 for x e [—2, 2] and p(x) = otherwise. 
(Hi) When c > -2, 

p(x) = i(6 2 x 2 + b )V a 2 — x 2 \x\ < a, 

7T 

= otherwise, 

where a 2 - ^^f- 2 ^ bo = c+ ^f+~\ and b 2 = \. 

When c e [0, oo), V is C 2 convex and V"(x) > 3 for \x\ > 1. In this case, The- 
orem 1.3 (hi) implies that W^MtiMv) ~^ with exponential convergence rate. When 
c < —2, [iy has two supports [—6, —a] and [a, b] which are disjoint. By Section 7.1 in 
Biane-Speicher [4], it is known that /j, t does not converge to [iy. See also Biane [3] and 
Remark 1.7. This also indicates that one cannot simultaneously prove a free version of 
the Holley-Stroock logarithmic Sobolev inequality and a free version of the Talagrand T 2 - 
transportation cost inequality under bounded perturbations of the distribution of eigenval- 
ues pn((Ix) — Zjj 1 n»<j \ x i ~ x j\ 2 UiLi e~ NV< - Xi) dx. Otherwise, by analogue of the proof of 
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Theorem 1.4 (iii), we may prove that pt converges to fiy with respect the W2-Wasserstein 
distance and hence in the weak convergence topology on ^(R). See also [25, 26] for a 
discussion on non-convex potentials. 

In the case c G [—2,0), as the global minimizer py of £y has a unique support, and 
all stationary point of py must satisfy the Euler-Lagrange equation Hp = \V ■, one can 
see that the Voiculescu free entropy £y has a unique minimizer py. As p t is the gradient 
flow of £y on ^ 2 (R), and since |S y (/i() = -2 J R [V'(x) - 2(Hp)(x)} 2 p(x)dx, we see that 
Ey(/j ( ) is strictly decreasing in time t unless pt achieves the (unique) minimizer py. This 
yields that Sy(/i t ) converges to some value. The question whether Wi{pt, f-iv) — > or even 
fit weakly converges to py as t — > oo for the above double-well potential V remains open. 
We would like to raise the following conjectures. 

Conjecture 6.1 Consider the double-well potential V(x) = jx 4 + |x 2 with c e [—2,0). 
Then p t converges to py with respect the W 2 -W asset stein distance and hence in the weak 
convergence topology on ^(R). 

In general, we may raise the following conjectures. 

Conjecture 6.2 Suppose that the potential V is a C 2 potential function with V"(x) > K\ 
for all \x\ > r and V"(x) > —K 2 for all \x\ < r, where K\,K 2 ,r > are some constants. 
Suppose further that £y has a unique minimizer which has a unique compact support. Then 
p(t) converges to py with respect the Wi-Wasserstein distance and in the weak convergence 
topology on ^(R). 

Conjecture 6.3 Suppose that the potential V is a C 2 potential function with V"(x) > 
—K for all x € R, where K > is a constant. Suppose further that Sy has a unique 
minimizer which has a unique compact support. Then p{t) converges to py with respect the 
Wi-Wasserstein distance and in the weak convergence topology on ^(R). 

Finally, let us mention the following conjecture due to Biane and Speicher [4]. 

Conjecture 6.4 Consider the double-well potential given by V(x) = ^x 2 + jX A , where g 
is a negative constant but very close to zero. Then p t weakly converges to py. 
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